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Abstract 

In this work our previous fully differential ionization cross section calculations using the semiclassical, impact param- 
eter method are improved by a new method suitable to calculate impact parameter values corresponding to different 
momentum transfers. This goal is achieved by two successive steps. First, using the transverse momentum balance dif- 
ferent projectile scattering angles are calculated for the binary and recoil peak regions as a function of the transferred 
momentum. Then, by treating the projectile scattering as a classical potential scattering problem, impact parameters 
are assigned to these scattering angles. The new method, which no longer contains empirical considerations, is tested 
calculating by fully differential ionization cross sections for single ionization of helium produced by fast C®^ ions. 
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In the last 5 years there was an intensive effort 
to measure fuUy differential cross sections by kine- 
matically complete experiments. First, the group of 
Schulz reported interesting data for the complete 
electron emission pattern in single ionization of 
helium by the impact of fast C^^ ions for certain 
momentum transfers [1,2]. The three-dimensional 
images were obtained using experimentally mea- 
sured fully differential cross section values on 
a cold-target-recoil-ion-momentum spectrometer 
(COLTRIMS) apparatus. Their results show a char- 
acteristic double-lobe structure with a binary peak 
and a smaller recoil peak. 

On theoretical side several calculations exist being 
able to reproduce the experimental data in the scat- 
tering plane (determined by the momentum of the 
scattered projectile and the momentum transfer vec- 
tors) [3,4]. At the same time, in the case of the per- 
pendicular plane there is no agreement between the 
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theoretical predictions and the experiment. Some 
authors have suggested that it may be important 
to include the internuclear interaction into the cal- 
culations [5]. On the other hand, very recently, the 
importance of taking into account the uncertainties 
of the experimental measurements and performing 
a convolution of the theoretical results on the exper- 
imental resolution [6,7] was proved. 

Previously, based on the semiclassical impact pa- 
rameter method we constructed a theoretical model 
to calculate fully differential cross sections for sin- 
gle ionization of helium by fast C^"*" ion impact [8] . 
Based on the semi-empirical version of the model a 
good agreement with the experiment was achieved 
in the scattering plane, while in the perpendicular 
plane a structure similar to that observed experi- 
mentally was obtained. It's main drawback was that 
the impact parameter values were selected based on 
the experimental data available in scattering plane. 

In this work the model is improved by the intro- 
duction of a new method for calculating impact pa- 
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rameter values corresponding to different momen- 
tum transfers. Consequently, the empirical part of 
the description may be eliminated and the applica- 
bility of the semiclassical impact parameter method 
to calculate fully differential cross sections may be 
confirmed. Then, the method is tested in case of sin- 
gle ionization of helium produced by ion pro- 
jectile with an energy of Eq ~ 100 McV/u. 

As described in [8] the semiclassical impact pa- 
rameter method may be used to calculate fully dif- 
ferential ionization cross sections. In this approxi- 
mation, the projectile is treated separately and it 
moves along classical trajectory. This implies that 
only the electron system needs to be; described by 
a time-dependent Schrodinger equation, while the 
projectile follows classical laws throughout. 

The ionization probability amplitudes for ioniza- 
tion of helium produced by fast C^^ ions are calcu- 
lated using first order time dependent perturbation 
theory. The initial state of the dielectronic system is 
described by a Hartree-Fock wavefunction [9] , while 
the final state of the system is described by symmet- 
ric combination of a hydrogenic and a continuum ra- 
dial wavefunction (calculated in the mean field of the 
final He"*" ion). Therefore, the ionization probabil- 
ity amplitude depending on the momentum transfer 
vector, ejected electron energy and electron ejection 
angles is reduced to a one-electron amplitude 
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where i and / represent the target system's initial 
and final electronic states, while the indices b and 
c represent bound and continuum states. Similarly, 
Ei and Ef are the energies of the corresponding (un- 
perturbed) states of the system and V denotes the 
time-dependent interaction between projectile and 
active electron. The projectile velocity is denoted by 
V and the integral is calculated through its classical 
trajectory along z axis. 

This amplitude is calculated expanding the final 
continuum-state wavefunctions into partial waves. 
As a result, amplitudes for transitions to ionized 
states with different angular momenta (a|^]„^) are 
obtained. 

The fully differential cross sections relative to the 
momentum transfer value q, ejected electron energy 
E and electron ejection angles 6 and ^ are obtained 
by the relation 
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where B is the impact parameter vector and I / and 
m/ are quantum numbers of the partial waves de- 
scribing the ejected electron. 

Certainly, the most sensitive part of the model 
is to assign impact parameter values to a certain 
momentum transfer. This task may be completed in 
two successive steps. 

First, the projectile scattering angle will be cal- 
culated by the use of the transverse momentum bal- 
ance [10] which states that the momentum transfer 
q is the sum of the transverse components of elec- 
tron's and residual ion's momenta. This vector rela- 
tion may be written in scalar form as 



Pti. = pI± + <f - ^Pe±q cos </) , 
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where pt_l is the transverse momentum taken by the 
residual ion and Pe± = Pe sin 9 is the transverse mo- 
mentum of the ionized electron. Further we assume, 
that the impact parameter is related to the momen- 
tum transfer to the residual ion, and take into ac- 
count the projectile-electron interaction separately. 

Furthermore, let us now suppose that the momen- 
tum transfer pT± is modifying only the direction 
of the projectile momentum vector. Supposing that 
we have to deal with small scattering angles (valid 
for rapid collisions) and neglecting the terms being 
many times smaller than 1 the projectile scattering 
angle may be calculated as 
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where po denotes the modulus of the projectile mo- 
mentum. 

In case of binary peak, one has to deal with = 

while in case of recoil peak the value of this angle is 
180°. In both cases the maxima are a.t 9 = 90°. For 
these peaks the scattering angles are 
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For the experimentally studied q = 0.75 a.u. and 
E = Q.5eV case these angles have values of 9^1^^^^ = 
4.2073 • 10-« rad and 9'^^lf = 1.0337 • 10"^ rad. 

Second, our next task is to assign impact param- 
eter values to the projectile scattering angles. Ac- 
cordingly, the projectile scattering will be treated as 
a classical potential scattering problem in the field 
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Fig. 1. Calculated projectile scattering angles for different 
impact parameter values. 

of the target helium system [11]. Therefore, the scat- 
tering angle may be calculated as 

where tq is the distance of closest approach defined 
by ^1 — — = and U{r) is the scattering 
potential. The simplest way to include the effect of 
the electrons around the target nucleus is to con- 
sider the potential to be a product of the Coulomb 
potential and the Bohr- type screening function [12] 

U{r) = ^P£2J^!^e-^ . (7) 
r 

If, particularly, the Thomas-Fermi potential [13] is 
used, then for the C^"*" -|- He colliding system the 
value of parameter a will be given by 

0.8853 

« - -vQ.23 (8) 
target 

expressed in a.u. 

Using this potential the integral (6) may be cal- 
culated numerically. Figure 1 shows the results of 
our calculations, namely the dependence of the scat- 
tering angle on impact parameter. From this data 
it turns out that impact parameters corresponding 
to previously calculated binary and recoil scattering 
angles are 2.47 and 0.68 a.u., respectively. 

Just for visualization purposes, based on our new 
theoretical results we have constructed the 3D im- 
age of the electron emission pattern for the studied 
case. The results presented on Figure 2 are calcu- 
lated using the previously determined impact pa- 
rameter values. The transition between these impact 
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Fig. 2. 3D imago of the electron emission pattern for sin- 
gle ionization of helium produced by 100 MeV/u C®+ pro- 
jectile impact calculated by semiclassical impact parameter 
method. The ejected electron energy is Ee = 6.5 eV and the 
momentum transfer is g = 0.75 a.u. 

parameters is realized smoothly in the < < 50° 
and 130° < 9 < 180°, 90° < (/) < 270° transition 
regions. As sketched on the figure, the initial projec- 
tile direction is along the z axis and the momentum 
transfer vector q is pointing nearly in x direction. 
In the 3D graph one can observe the presence of the 
characteristic double-lobe structure towards the x 
axis with binary peak at 6* = 90°, = 0° and recoil 
peak at 6* = 90°, (j) = 180°. The electron emission 
pattern is almost identical with our previous result 
obtained by the semi-empirical model [8] . 

In order to analyze in detail the obtained results, 
cross section values for scattering plane and perpen- 
dicular plane are plotted separately in figure 3. Thus, 
the results of the present theory may be compared 
in absolute value to experimental data, previous re- 
sults of the semi-empirical version of the model and 
CDW calculations. 

The top panel of the figure shows fully differential 
cross sections in the scattering plane for different 
electron ejection angles <j). The present model gives 
good agreement with experimental data of Schulz et 
al. [1]. However, a small difference may be detected 
between the current version and the semi-empirical 
version of the model mainly in the binary peak re- 
gion. The difference may be explained knowing that 
in case of the previous semi-empirical model impact 
parameters of 2.2 and 0.7 a.u. have been selected 
based on the available experimental data [8] . Even 
so, the difference between results looks to be smaller 
than the experimental errors. 
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Fig. 3. Theoretical results in scattering (top) and perpendic- 
ular (bottom) planes compared with experiments [6] for the 

same case as in figure 2. The solid curve shows the present 
theory, dashed line represent results of the semi-empirical 
version of the model [8], while the dash-dotted line is ob- 
tained by the CDW model [6]. 



decessor, the model describes well the fully differ- 
ential cross sections for relatively small momentum 
transfer values. The characteristic structures in the 
perpendicular plane have also been reproduced, dis- 
crepancies with the experiments are only in the mag- 
nitude of the cross sections, which may be explained 
by the experimental uncertainties. 

Finally, it has to be mentioned, that the main goal 
of this work was not to improve the previous results 
of the scmiclassical impact parameter model used 
to calculate fully differential cross sections, but to 
create a pmc theoretical background for it. Accord- 
ingly, all previous results of the model have been 
confirmed, and we have given a theoretical method 
for obtaining the impact parameter used in the cal- 
culations. 

The present work has been supported by the Ro- 
manian National Plan for Research (PN II), contract 
No. ID 539. 



In case of the perpendicular plane almost the same 
results as the semi- empirical ones have been ob- 
tained. The bottom panel of figure 3 shows fully 
differential cross sections in this plane for different 
electron ejection angles 9. The curve shows the same 
behavior as the experimental data with strong max- 
ima at 6* = 80° and = 280°. A third smaller maxi- 
mum is also obtained at direction of 9 = 180°. 

Here it has to be noted that, compared to the 
isotropic results of the CDW model, a better agree- 
ment in shape has been obtained by the present 
metod. However, the magnitude of the cross section 
is smaller than the experimental one. The recently 
reported inclusion of the experimental momentum 
uncertainties [6] should also improve the agreement 
between theory and experiments by increasing the 
cross sections in the perpendicular plane. Our result 
is consistent with the conclusions of Diirr et al. [7] , 
that the experimental uncertainties are responsible 
only partly for the structure observed in the perpen- 
dicular plane, half of the value of the maxima may 
be due to some real physical effect. 

In conclusion, the theoretical model based on the 
first order, semiclassical, impact parameter approxi- 
mation for describing kinematically complete exper- 
iments has been revised by eliminating any empir- 
ical considerations from it. This has been achieved 
by calculating impact parameter values correspond- 
ing to different momentum transfers. Then, the new 
model was applied for studying single ionization of 
helium by impact with fast ions. As its pre- 
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